J/ip suppression and p T spectra in RHIC and LHC energy collisions 
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In a hydrodynamic model, we have studied J/ip production in Au+Au/Cu+Cu collisions at 
RHIC energy y / s=200 GeV. At the initial time, J/ip's are randomly distributed in the fluid. As the 
fluid evolve in time, the free streaming J/ip's are dissolved if the local fluid temperature exceeds 
a threshold temperature Tj/^,. Sequential melting of charmonium states (x c , ip 1 and J/ip), with 
melting temperatures T Xa = T^i ~ 1.2T C , Tj/^ ~ 2T C and feed-down fraction F « 0.3, explains 
the PHENIX data on the centrality dependence of J/ip suppression in Au+Au collisions. J/ip pr 
spectra and the nuclear modification factor in Au+Au collisions are also well explained in the model. 
The model however over predict centrality dependence of J/ip suppression in Cu+Cu collisions by 
20-30%. The J/ip pr spectra are under predicted by 20-30%. The model predict that in central 
Pb+Pb collisions at LHC energy, -^5=5500 GeV, J/ip's are suppressed by a factor of ~ 10. The 
model predicted J/ip pr distribution in Pb+Pb collisions at LHC is similar to that in Au+Au 
collisions at RHIC. 

PACS numbers: PACS numbers: 25.75.-q, 25.75.Dw 



I. INTRODUCTION 

High energy heavy ion collisions give the opportunity 
to study QCD matter at high density and temperature. 
Recent experiments at Relativistic Heavy Ion Collider 
(RHIC) at Brook Haven National Laboratory, give strong 
indications that in Au+Au collisions, a collective QCD 
matter is formed [l|, [2j, [3J, Ij] • However, whether the mat- 
ter can be qualified as the lattice QCD predicted Quark- 
Gluon Plasma Q is still a matter of debate. J/ip sup- 
pression is long recognized as an important tool to iden- 
tify the possible phase transition to quark-gluon plasma. 
Matsui and Satz [f| , predicted that in presence of quark- 
gluon plasma (QGP), binding of a cc pair into a J/ip 
meson will be hindered, leading to the so called J/ip sup- 
pression in heavy ion collisions [6(. It was later realised 
that J/ip's are not dissolved just at the critical tempera- 
ture T c . Lattice based potential models indicate that dis- 
sociation temperature for J/ip is Tjm « 2.1 T c |7|. The 
excited states \ c and ipi on the other hand are dissolved 
close to the critical temperature T x « 1.2T C , T^, « 1.1T C 
0. 

PHENIX collaboration have measured J/ip yield in 
p+p [1, J! [nj, d+Au [I, OH, Au+Au pj, [ll and 
Cu+Cu [14j collisions at RHIC. J/ip's are suppressed 
in d+Au, Au+Au and Cu+Cu collisions. Suppression 
in d+Au collisions is due to a number of cold nuclear 
matter (CNM) effects, e.g. gluon shadowing, gluon en- 
ergy loss, J/-0-nucleon inelastic collision. J/ip suppres- 
sion in hot nuclear matter is measured in Au+Au and 
Cu+Cu collisions. In Au+Au/Cu+Cu collisions, sup- 
pression increases with centrality of the collisions and in 
very peripheral collisions, suppression is consistent with 
CNM effects. Also J/ip's are more suppressed in for- 
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ward rapidity than at mid rapidity. Measurements also 
indicate that models [HI, [l6| , which predict for J j ip en- 
hancement in RHIC energy collisions, are not consistent 
with experiments. On the other hand, phenomenologi- 
cal threshold model [Tt], EH, which mimic the melting 
of J/ip's in a QGP medium, successfully explain the sup- 
pression in Au+Au collisions but not in Cu+Cu collisions 
[19L I20I |2~l| . The threshold model do not account for the 
possible recombination of cc pairs. In 22] PHENIX data 
were analysed in the comover model. It was claimed that 
the comover model , including recombination, is consis- 
tent with the centrality dependence of J/ip suppression 
both in Au+Au and Cu+Cu collisions [23]. However, we 
noticed that in mid-rapidity Au+Au collisions, the model 
do not reproduce the sudden increase of suppression be- 
yond N par t ~ 150. 

Recently, following Gunji et aL[23j, we have developed 
a 'Hydro+J/'0' model to simulate J/ip production in nu- 
clear collisions [24j . The QGP fluid evolves in 2+1 di- 
mensions. At the initial time, J/ip's are randomly dis- 
tributed in the transverse plane. As the fluid evolve, free 
streaming J/ijj's arc melted if the local fluid tempera- 
ture exceeds a critical value. A considerable percentage 
of observed J/ip's are from decay of higher states, \ an d 
tp/ [25] ]. The experimental J/ip suppression pattern in 
mid-rapidity Au+Au collisions is well explained by se- 
quential melting of Xc 4> f and J/ip in the dynamically 
expanding fluid. The estimated melting temperatures, 
Tj/ip — 2.1T C , T Xc = = 1.1T C are in agreement with 
the lattice motivated calculations [7|. The fraction of 
the higher states (x c + 4>') is estimated to be F = 0.3. 
The model however fails to reproduce the centrality de- 
pendence of J/ip suppression in Cu+Cu collisions. Sup- 
pression is over predicted. In the present paper, we have 
improved upon the model and studied transverse momen- 
tum distribution of J/ip's in Au+Au and in Cu+Cu colli- 
sions. The improved model, in addition to the centrality 
dependence of J / ip suppression, well reproduces the J / ip 
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Pt spectra in Au+Au collisions. However, in Cu+Cu 
collisions, the improved model over predict J /if) suppres- 
sion. The pt spectra is under predicted. We have also 
used the model to predict for the centrality dependence 
of J/ip suppression and pt spectra in Pb+Pb collisions at 
LHC energy. In the model, J ftp's are greatly suppressed 
in LHC energy. 

The plan of the paper is as follows: in section [Til we 
briefly describe the 'Hydro+J/^' model. In section HIT1 
we have analysed PHENIX data on the centrality depen- 
dence of J ftp suppression, J / ip px spectra in Au+Au and 
Cu+Cu collisions. In section ITVl we give predictions for 
centrality dependence of J / ip suppression in Pb+Pb colli- 
sions at LHC. Predictions for J ftp pr spectra, in 0-20%, 
20-40%, 40-60% and 60-92% Pb+Pb collisions at LHC 
energy are also given in section IIVI Lastly, summary 
and conclusions are drawn in section fVl 



II. THE HYDRODYNAMIC MODEL FOR J/ip 
SUPPRESSION 



Details of the hydrodynamic model used here can be 
found in 0, Briefly, it is assumed that in high 

energy nuclear collisions, a deconfined phase (QGP) is 
produced, which expands, cools, undergoes 1st order 
phase transition to hadronic fluid at the critical temper- 
ature (T c =164 MeV) and then further cools to freeze- 
out at temperature Tf=130 MeV. Assuming longitudi- 
nal boost-invariance, the space-time evolution of the fluid 
is obtained by solving the energy-momentum conserva- 
tion equation d^T^ — [26]. Assumption of boost- 
invariance limits the model to mid-rapidity region only. 
J ftp's are assumed to be produced in initial NN collisions. 
They are free streaming, unless dissolved in the QGP 
medium. Hydrodynamic models require energy density, 
fluid velocity distributions at the initial time Tj . The ini- 
tial energy density of the fluid in the transverse plane 
was parameterised as 0, HH , 



e(h, x, y) = e [0.75Ap art (b, x, y) + .25N coU (h, x, y)} (1) 

where N part (h,x,y) and N co u(h,x,y) are the transverse 
profile for the participant number and binary collisions 
number in an impact parameter b Au+Au collision. The 
initial fluid velocity was assumed to be zero, v x (h, x, y) — 
v y (h,x,y) = 0. The constant e depend only on the col- 
lision energy and not on centrality of the collisions. The 
initial time Ti and the constant Eq are chosen to repro- 
duce the pt distribution of identified particles in central 
Au+Au collisions. For b=0 Au+Au collisions, it corre- 
spond to central energy density, £=30 GeV/fm 3 or cen- 
tral entropy density Si n i — 110 /to -3 , at the initial time 
Ti=0.6 fm/c. Solution of hydrodynamic equations require 
an equation of state. We have used the equation of state 
EOS-Q [26j, incorporating 1st order phase transition with 
critical temperature T c =164 MeV. 



To obtain the survival probability of J ftp's in an ex- 
panding medium, we proceed as follows: at the initial 
time Ti — 0.6/m/c, we randomly distribute a fixed num- 
ber of J ftp's in the transverse plane. They are assumed 
to be free streaming [271 ] unless dissolved in the medium. 
Each J ftp is characterised by 4 random numbers. Two 
random numbers (i?i,i?2) indicate its transverse posi- 
tion (rjj, and two random numbers (i?3,i?4) its trans- 
verse momentum px- In [H, [24J, J ftp's are spatially 
distributed according to initial fluid temperature pro- 
file. However, initial NN collisions produces cc pairs, 
which later evolve to J ftp. The formation time of J ftp ~ 
0.5 fm, incidentally is similar to the thermalisation time 
of the QGP fluid. Then, spatial distribution of initial 
J ftp's may not follow the spatial distribution of the fluid 
temperature, rather it is more likely follow the trans- 
verse profile of the binary collision number. Presently, 
we distribute the random numbers R\ and i?2 according 
to the transverse profile of the number of binary collisions 
(N coll ). 

PHENIX collaboration have measured Jf ip pt spectra 
in p+p collisions at RHIC The power law [H, 
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with A — 4.23 and B = 4.1 rather well describe the 
invariant distribution of measured J ftp's in p+p collisions 
. The random number R 3 is distributed following the 
power law EqJ^l The random number is distributed 
uniformly within [0-27r]. 

The survival probability of a J ftp inside the expanding 
QGP is calculated as [23j |. 



S J/->p( T ) = ex P 



T dis (T(r ± (T')))dT' 



(3) 



where T(r±) is the temperature of the fluid at the trans- 
verse position r±, Tdi S (T) is the decay width of J ftp at 
temperature T. Ti is the initial time for hydrodynamic 
evolution. We continue the evolution till the freeze-out 
temperature Tp=130 MeV. In our earlier work [24J, for 
the decay width Tdis we had used the simplest prescrip- 
tion, 
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The decay width EqQ] is rather rough. At the melting 
temperatur e T / /,;,, it abruptly changes from zero to oo. 
Gunji et al (Hf, in their analysis PHENIX Au+Au data 
used the following form for the decay width, 



Fdis{T) — oo; T > Tj/^ 

T dls (T) = a(T/T c -l) 2 ; 
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(5) 
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FIG. 1: (color online) The filled circles are the PHENIX data 
[13 . l33l . l34l ] for the centrality dependence of nuclear modi- 
fication factor (Raa) for J /if) in mid- rapidity Au+Au colli- 
sions. 'Hydro+ J/ip' model predictions for four values of the 
decay width controlling parameter a=0.1, 0.2, 0.3, and 0.4 are 
shown in 4 panels. In each panel, the red and green lines are 
Raa for direct J ftp's and excited states ip/ \ c . The melting 
temperatures are T J/i} = 2T C , T^, = T Xc = 1.2T C . The black 
lines are Raa assuming 30% feed-down. 



In EqO a is the thermal width of the state at T/T c =2. 
NLO perturbative calculations suggest that a > QAGeV 
[29[. Gunji et al ^ observed that for a > 0.4, the 
PHENIX data are not well described. Data are best de- 
scribed with a=0.1. We note that Eqj5]do not take into 
account of regeneration of charmonium states. Lattice 
QCD calculations [3(| indicate D- and B-like states can 
exist in sQGP. D- and B-like states in sQGP can have 
direct impact on J/tp production, facilitating regenera- 
tion. They can also provide for the resonant cross sec- 
tion for heavy quarks and play an play an essential part in 
thermalisation and collective flow of charm and bottom 
quarks [3l| . 

Any model of J/tp suppression must account for the 
experimental observation that a substantial fraction of 
the measured J ftp's are from decay of the excited char- 
monium states Xc and ip' [111- To calculate the survival 
probability for the excited states \c an d ip', we use the 
same procedure as described above for the ground state 
J/tp- Above T Xa and T^ h the excited states Xc and ip' 
are assumed to melt. Further noting that the lattice mo- 
tivated calculations indicate T Xc « Tm, we define a 
common temperature T x = T Xc — T^,, above which all 
the states Xc and ip' are dissolved. For feed-down fraction 
F, the J/tp survival probability is then obtained as, 



S QGP = (l-F)Sj^+FS x (6) 
As mentioned earlier, PHENIX collaboration in d+Au 
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FIG. 2: (color online) Black and red circles are PHENIX 
data for pr dependence of nuclear modification factor (Raa) 
for J/tp's in 0-20%, 20-40% 2 40-60% and 60-90% centrality 
Au+Au [JJj] and Cu+Cu [l4| collisions respectively. Er- 
rors are statistical only. The black and red lines are the 
'Hydro+ J/tp' model predictions for Au+Au and Cu+Cu col- 
lisions. Decay width controlling parameter a=0.2. CNM ef- 
fects are included (see text). 



collisions has studied cold nuclear matter effect on J/tp 
suppression [ij J ftp's are suppressed in d+Au col- 
lisions also. The suppression is consistent with Glauber 
model of nuclear absorption with J/^-nucleon absorption 
cross-section <r a b s — 2 ± 1 mb [H[ . If cold nuclear matter 
effect is taken into account, the survival probability of 
J ftp can be obtained as, 



5' 
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S, 



QGP 



X 5, 



CNM, 



(7) 



where Scnm is the survival probability in cold nuclear 
matter. 



III. J/tp SUPPRESSION IN RHIC ENERGY 

A. Au+Au collisions 

Experimental survival probability for J/ijj in mid- 
rapidity Au+Au collisions are shown in Fig[T] where 
gluon shadowing and nuclear absorption with cr a bs=l mb 
is taken into account as cold nuclear matter (CNM) ef- 
fect [13, [H, [IH • In FiglU we have shown the results of 
the present 'Hydro+ J/ip' model. They are obtained for 
melting temperatures Tjiy = 2T Cl T^, — T Xa — 1.2T C . 
Results are shown for four values of the decay width con- 
trolling parameter a=0.1, 0.2, 0.3 and 0.4. The red and 
green lines are survival probability for the states J/ip and 
X- The black line is the combined suppression with 30% 
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FIG. 3: (color online) The colored circles are the PHENIX 
data H3 for J/ip p T distribution in 0-20%, 20-40%, 40-60% 
and 60-90% centrality Au+Au collisions. The black lines are 
the J/ip pt spectra in , J/ip+ hydro' model predictions. The 
CNM effect is included (see text). 
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FIG. 4: The black circles are PHENIX data [3 on the cen- 
trality dependence of J/ip suppression in Cu+Cu collisions. 
The black line is the 'Hydro+ J/ip' model prediction for the 
J/ip suppression in Cu+Cu collisions. 



feed-down. Sequential melting of charmonium states ex- 
plains the PHENIX data on the centrality dependence 
of J/ip suppression in Au+Au collisions. As in [23j |. 
J/ip suppression depend on the decay width controlling 
parameter a, increasing with a. However, the depen- 
dence seems to be smaller than found by Gunji et al [2a |. 
For example, in the most central collision (N part = 325) 
changing a from 0.1 to 0.4 changes Raa by less than a 
few percent. Good description of the data for cv=0. 1-0.4, 
leave little or no room for charmonium regeneration. 

PHENIX collaboration have measured J/ip transverse 
momentum distribution and also the nuclear modifica- 
tion factor (R AA (p T )), in 0-20%, 20-40% , 40-60% and 
60-92% centrality Au+Au collisions In FigEJ in 

four panels, PHENIX results for the nuclear modifica- 
tion factor (Raa) in 0-20%, 20-40%, 40-60% and 60-92% 
centrality Au+Au collisions are shown. Suppression in- 
creases as the collision centrality increases. Also within 
the errors, suppression is approximately px independent. 
Note that PHENIX data on Raa(pt) are not corrected 
for the CNM effects. We have approximated CNM ef- 
fect by the Glauber model of nuclear absorption with 
32|. In FiglH black lines are 'Hydro+J/^' 
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model prediction. The decay width controlling parame- 
ter is ev=0.2. Sequential melting of charmonium states 
correctly reproduces the PHENIX data on the px depen- 
dence of nuclear modification factor. Only in very pe- 
ripheral (60-92% centrality) collisions, the model do not 
reproduce the experimental data in px range 1-3 GeV. 
It is not surprising. Hydrodynamics models are not very 
accurate at large pt or in very peripheral collisions [26j . 
In FigjJH we have compared model predictions against the 
PHENIX data on J/tp p T distribution in 0-20%, 20-40%, 
40-60% and 60-92% centrality Au+Au collisions. As ex- 



pected, the 'Hydro+ J /ip J model, with sequential melting 
of charmonium states, also explain the PHENIX data on 
the J 1 ip Pt spectra. 

PHENIX collaboration also measured the mean square 
transverse momentum (< p T >) for J/ip. In table HI 
PHENIX measurements of < p T > in different centrality 
ranges of Au+Au collisions are shown. Within the errors, 
< p\ > do not show variation with collisions centrality 
and is approximately same as in pp collisions. As shown 
in the table, 'Hydro+ J/tp' model also do not show any 
appreciable variation of < p\ > with collision centrality. 
FigsQ] [H [3] and table Q] clearly indicate that sequential 
melting of charmonium states, with melting temperature 
Tj/,0 = 2T C , T^), = T Xa = 1.2T C and feed down fraction 
F=0.3, explains most of the features of J /ip production 
in Au+Au collisions at RHIC energy. 



TABLE I: PHENIX results for J/tp < p T > are compared 
with 'Hydro+ J/ip' model calculations. Melting temperatures 



of J /if) and x are Tj/^ 
fraction is F=0.3. 



2T C and T Y = 1.2T C . The feed down 
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FIG. 5: (color online) The colored symbols are the PHENIX 
data for J/ip p T distribution in 0-20%, 20-40%, 40-60% and 
60-90% centrality Cu+Cu collisions [14|. The black lines are 
the J/ip pt spectra in 'Hydro+ J/ip'' model. The melting tem- 
perature are Tj/j, = 2.0T C and T x = T Xc = 7> = 1.2T C . The 
feed-down fraction is F = 0.3. 
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FIG. 6: 'Hydro+ J/ip' model predictions for the centrality 
dependence of J/ip suppression in Pb+Pb collisions at LHC 
energy. The solid line is the model prediction with the cold 
nuclear matter effects. For comparison, PHENIX data on the 
centrality dependence of J/ip suppression in Au+Au collisions 
13, 33, 34] arc also shown. 



B. Cu+Cu collisions 



Recently PHENIX collaboration published their anal- 
ysis of J/ip measurements in Cu+Cu collisions [HI]. In 
FigHl PHENIX data [3 on the centrality dependence 
of J/ip suppression are shown (the black circles). Data 
are not corrected for the cold nuclear matter (CNM) ef- 
fects. In the most central collision, J/ip's are suppressed 
by a factor of ~ 2. Survival probability of J/ip increases 
as the collision centrality decreases till N part =lQ.7 . In 
more peripheral (N part =6A), collisions, survival proba- 
bility decreases again. The decrease in survival probabil- 
ity as participant number decreases from N part —16.7 to 
6.4 is interesting. All the theoretical models predict con- 
tinuous increase of suppression as the collision centrality 
increases. 

The black line in Fig[J]is the 'Hydro+J/?/'' model pre- 
diction for the centrality dependence of J/ip suppres- 
sion in Cu+Cu collisions. The model parameters are 
same as it was for Au+Au collisions, i.e. Tji^ = 2T C , 
T$i = T Xc = 1.2T C , F = 0.3 and a=0.1. We have approx- 
imated the CNM effects by the Glauber model of nuclear 
absorption with J/f/;-nucleon absorption cross section 
a a ft s =2 mb. Sequential melting of charmonium states 
fails to explain the PHENIX data on the centrality depen- 
dence of J/ip suppression. In central collisions, the model 
over predict the suppression by 20-30%. As discussed ear- 
lier, we have neglected charmonium regeneration. If re- 
generation of J/ip is important in Cu+Cu collisions, the 
gap between theory and experiment may be filled by the 
regenerated charmoniums. However, as shown earlier, in 
Au+Au collisions, there is little or no scope for recombi- 
nation of cc pairs in 'Hydro+J/^' model. It is unlikely 



that regeneration of charmoniums will be important in a 
smaller system like Cu+Cu. Indeed, in [22j], where J/ip 
suppression in Au+Au and Cu+Cu collisions were anal- 
syed in the comover model, including recombination, it 
was observed that recombination effect reduces suppres- 
sion by ~50% in central Au+Au collisions, but reduction 
is only ~ 10% in central Cu+Cu collisions. 

PHENIX collaboration also measured J / ip pr distribu- 
tion in Cu+Cu collisions [lj|. In FigO we have shown 
the PHENIX measurements for J/ip px spectra in 0- 
20%, 20-40%, 40-60% and 60-92% centrality Cu+Cu col- 
lisions. In Figf51 the black lines are the pr spectra in 
'Hydro+ J/ip' model with sequential melting of charmo- 
nium states. CNM effects are included. PHENIX data 
are reproduced within 20-30%. Quality of fit to the data 
is better judged if nuclear modification factor Raa(j>t) 
for J/ip's are compared . In Fig[2l the red symbols are 
Raa(pt) for J/ip in Cu+Cu collisions (3f| and red lines 
are the model predictions. Compared to Au+Au colli- 
sions, J/ip's are slightly less suppressed in Cu+Cu colli- 
sions. Also, as it is for Au+Au collisions, in Cu+Cu col- 
lisions also, Raa(pt) is approximately px independent 
and increases with centrality . In 0-20%, 20-40% and 40- 
60% centrality Cu+Cu collisions, sequential melting of 
charmonium states over predict the data by 20-30%. Ev- 
idently, charmonium production mechanism in Cu+Cu 
collisions is at variance with the charmonium production 
in Au+Au collisions. 
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FIG. 7: Black lines are the 'Hydro+ J/ip' model predictions 
for the J/ip p T distribution in 0-20%, 20-40%, 40-60% and 
60-92% centrality Pb+Pb collisions at LHC. For comparison, 
PHENIX data for Au+Au collisions are also shown. 20- 
40% centrality collision onwards, the curves (and data) are 
progressively divided by a factor of 10. 



IV. J/ip SUPPRESSION IN PB+PB COLLISION 
AT LHC 

'Hydro+ J/ip 1 model can be used to predict for the 
J/ip suppression in Pb+Pb collisions at LHC energy 
(-^=5500 GeV). The initial conditions for hydrodynamic 
evolution in Pb+Pb collisions at LHC energy is unknown. 
We have used the existing knowledge at RHIC energy to 
extrapolate to LHC energy and obtained the initial con- 
dition for the hydrodynamic evolution. Central tempera- 
ture of the fluid in b=0 Pb+Pb collisions is 7^=421 MeV 
(for details see 24]). This can be contrasted with cen- 
tral temperature Ti=351 MeV in b=0 Au+Au collisions 
at RHIC. The 'Hydro+ J/ip 1 model also require J/ip pr 
distribution in p+p collisions as an input. It is yet to be 
measured. While there are many theoretical predictions 
for J/ip pt distribution at LHC energy, most of them 
are in the framework of pQCD and predicts for large pr- 
'Hydro+J/^' model on the other hand is suitable for low 
Pt particles. Presently, we approximate the J/ip px dis- 
tribution in p+p collisions at LHC energy by the same 
power law as it is for RHIC collisions (Eq(2J). Particle 
multiplicity increases from RHIC to LHC energy. As- 
suming logarithmic increase with energy, at LHC energy, 
J/ip multiplicity can be increased by a factor of ^1.6 [24| . 
To be consistent with the increased multiplicity at LHC, 
we normalise Eq[2]by the factor of 1.6. 

At LHC, initial fluid temperature is high (T;=421 
MeV) and thermal production of charm may not be neg- 
ligible. Indeed, simple estimate indicate that from RHIC 
to LHC, thermal charm production can increase by a fac- 
tor of ~ 2. Recombination of thermally produced charm 
and anti-charm quarks at LHC may not be negligible. 



Indeed, statistical hadronisation model [3(| favors strong 
recombination effect in central Pb+Pb collisions and pre- 
dict enhancement of J/ip rather than suppression. How- 
ever, we continue to neglect the recombination effect at 
LHC energy. As shown earlier, in the Hydro+ J/ip model, 
existing Au+Au data do not favor recombination and the 
parameters controlling recombination effect at LHC en- 
ergy can not be determined. The predictions for J/ip 
suppression given below then can be considered as an 
upper bound for suppression at LHC. 

In Figj6l 'Hydro+ J/ip' model prediction for the cen- 
trality dependence of J/ip suppression in Pb+Pb colli- 
sions is shown. The melting temperature of the states 
and the feed-down fraction are unchanged, Tj/^ = 2T C 
and T x = 1.2T C and F=0.3. They are property of the 
deconfined medium and do not depend on the collision 
energy. In Fig[|5J the solid line is 'Hydro+J/^' model 
prediction for the survival probability in Pb+Pb colli- 
sions at LHC. We have neglected the CNM effect here. If 
CNM effects are included, the suppression will increase 
further. For comparison, PHENIX data for Au+Au col- 
lisions at RHIC are also shown in Fig[6j Except for very 
peripheral collisions, J/ip are more suppressed at LHC 
energy than at RHIC energy. In central collisions, J/ip's 
are suppressed by a factor of ^10. At LHC energy, the 
fluid temperature is high and melting conditions for char- 
monium states are easily met than in Au+Au collisions 
at RHIC. Large suppression in LHC energy collisions will 
have experimental consequences. As stated earlier, from 
RHIC to LHC, particle multiplicities are expected to in- 
crease by a factor of ~ 1.6. The J/ip multiplicity is also 
expected to increase by a similar factor. But in central 
collisions from RHIC to LHC energy, J/ip suppression 
increases by a factor ~ 2.5. Then effectively J/ip yield 
will be reduced in LHC collisions by a factor of ~ 1.5. 
For similar statistics as in RHIC Au+Au collisions, ex- 
periments have to run for longer duration at LHC. 

Lastly, in FigJTJ we have shown the simulated J / ip pt 
spectra in 0-20%, 20-40%, 40-60% and 60-92% centrality 
Pb+Pb collisions. For comparison, we have also shown 
the PHENIX data for J/ip pt spectra in Au+Au colli- 
sions dHHH- It appears that J/ip pr spectra in Pb+Pb 
collisions at LHC is similar to that in Au+Au collisions 
at RHIC, may be slightly reduced in central collisions. 



V. SUMMARY AND CONCLUSIONS 

To summarise, in an improved version of 'Hydro+ J/ip' 
model [24j , we have studied J/ip suppression and px 
spectra in mid-rapidity Au+Au and Cu+Cu collisions 
at RHIC. It is assumed that in the collisions, a decon- 
fined phase (QGP) is produced, which expands, cools, 
undergoes 1st order phase transition to hadronic fluid at 
the critical temperature (T c ) and then further cools to 
freeze-out. The space-time evolution of the fluid is ob- 
tained by solving the hydrodynamic equations for ideal 
fluid in 2+1 dimensions [26| . J/ip's are assumed to be 
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produced in initial NN collisions. They are randomly 
distributed according to transverse profile of the binary 
collision number. As the fluid evolve in time, the "free 
streaming" J/ip's are melted if the local fluid tempera- 
ture exceed a critical temperature Tj/^. Similarly, the 
states Xc and ipf are assumed to melt above a critical 
temperature T x = T^p, = T Xc . Sequential melting of 
the charmonium states, ipt, x an d J/ip, above critical 
temperatures: Tji^ = 2T C , T^, = T Xc = 1.2T C with feed- 
down fraction F = 0.3, well explain the PHENIX data on 
the centrality dependence of J/ip suppression in Au+Au 
collisions. PHENIX data on J/ip pt spectra in 0-20%, 
20-40%, 40-60% and 60-92% centrality Au+Au collisions 
are also well explained. Analysis suggests that sequen- 
tial melting of charmonium states, in an expanding QGP 
fluid, is consistent with the PHENIX data on J/ip sup- 



pression and J/ip pt distribution in Au+Au collisions at 
RHIC. Sequential melting of charmonium states on the 
other hand appears to be inconsistent with experimental 
data on J/ip production in Cu+Cu collisions. Central- 
ity dependence of J ftp suppression is over predicted by 
20-30%. The J/ip pr spectra in Cu+Cu collisions are 
under predicted also by 20-30%. Within the model ap- 
proximation, J/ip production mechanism in Cu+Cu col- 
lisions is at variance with the production mechanism in 
Au+Au collisions. We have also predicted for the central- 
ity dependence of J/ip suppression and J/ip pt spectra 
in Pb+Pb collisions at LHC. The model predicts that in 
central Pb+Pb collisions, J/tp's are suppressed by a fac- 
tor of ~ 10. The J/ip pr spectra appears to be similar 
to that in Au+Au collisions at RHIC. 
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